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Abstract It is anticipated from hierarchical clustering theory that there are scal- 
ing relationships between halos over a wide range of mass. Observationally it can 
be difficult to identify the markers that characterize these relationships because of 
the small numbers of visible probes and confusion from contaminants in projec- 
tion. Nonetheless, in favorable circumstances it is possible to identify a very useful 
marker: the radius of the caustic at second turnaround. In a few favorable circum- 
stances it is possible to identify the radius of first turnaround, or zero velocity sur- 
face about a collapsed region. It will be shown that specifically the radius of second 
turnaround scales as anticipated over three orders of magnitude in mass from 10 12 to 
10 15 M . Halos are characterized by zones of dispersed velocities within the second 
turnaround caustic and zones of infall between the first and second turnaround radii. 
The inner zone is populated in the majority by gas poor morphologies and the outer 
zone is populated in the majority by gas rich morphologies. The numbers of dwarfs 
within the inner zone is roughly constant per unit halo mass. 



1 Introduction 

In the imaginary world of simulations, researchers have a well developed picture 
of the collapse of matter into halos. Over time, small halos are absorbed into larger 
units. Collapsed regions filled with substructure can be defined by hundreds or thou- 
sands of particles. Halos can be identified with precision within those simulations 



(Bovlan-Kolchin et al., 2009) 



In the real world, most galaxies are observed to lie in groups or clusters. How- 
ever, membership may be so limited that the structure is ill defined. There is no 
concensus among observers about what is meant by the terms 'group' and 'clus- 
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ter'. For example, it is commonly accepted that we live in something called the 
Local Group. A quantitativ e boundary that would inc lude the traditional members 



is the zero-velocity surface dKarachentsev et al. ■ 120091) . Galaxies inside this surface 



are infalling or on more complex bound orbits. Galaxies outside this surface are 
participating in the Hubble expansion. But now consider the Virgo Cluster. This 
entity would traditionally be taken to include the regio n within 2 Mpc where sev- 
eral thousand galaxies are following randomized orbits (IBinggeli et al.L 1 19871) . The 
zero-velocity surface around the Virgo Cluster lies at a radius of ~ 7 Mpc, almost 
half of our distance of 16.7 Mpc from the cluster. The implicit traditional definitions 
of the Local Group and the Virgo Cluster are inconsistent. 

It does not seem useful to draw a distinction between the terms 'cluster' and 
'group'. Clusters contain a lot of galaxies and groups contain only a few, but there 
is no clear demarcation between the two. In this article, the two terms will be used 
interchangeably. 

So how should an observer define a group? Structure forms as sufficient mat- 
ter accumulates through gravitational attraction to decouple from the expansion of 
the universe. The matter falls together on nearly radial orbits and, given enough 
time, evolves toward dynamic equilibrium. Two natural dimensions to describe a 
collapsed structure are the gravitational or virial radius, r g and the radius at 200 
times the 'critical' density for a matter-dominated closed universe, r2oo. Here, lower 
case r implies a 3-dimensional radius and an upper case R will be used to denote 
a projected radius. Accordingly, if all galaxies in a group sample are given equal 
weight (ie, galaxies of whatever luminosity are considered to be test particles within 
an environment dominated by distributed dark matter) then the virial radius is de- 
fined as 

N 2 
Li<jt/Rij 

where Rjj is the projected distance between pairs, and a virial mass estimate is given 
for the group by 

M v = atfr g /G = (n/2)a<7pR g /G (2) 
where the line-of-sight velocity dispersion for the sample of N galaxies is 



J>- < v >) 2 /N . (3) 

The group mean velocity is < v > and a — 3 if orbits are isotropically distributed. 

The alternative dimension r2oo is defined to coincide with a density 200p CT „ 
where p cnr = ml/^TzG atredshift zero. This parameter can be determined in simu- 
lations with large numbers of test particles but it is not such a useful con struct in the 
context of observations of small groups. However. iMahdavi et al. d 19991) calculate a 



ratio between r2oo and r g assuming M(r) <*; r that can be reduced to the relationship 
(atz = 0) 
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There may be other dimensions that are observationally useful. Consider the 
collapse of a spherically symmetric ove rdense region in the expanding universe 
dGunn & Gotdll972tlBertschingerlll985l) . The time of collapse, t c , depends on den- 
sity, p: t c o= p -1 / 2 . in the approximation of spherical collapse, all across the universe 
at t c = today overdense regions with a common density will be separating from the 
cosmic expansion, creating zero-velocity surfaces. The radius that encloses one of 
these regions will be called the first turnaround radius, r\ t . 

At earlier times, regions of successively higher density collapsed. There was a 
time and a corresponding density that lead to collapse and re -expansion to a pause 
befor e recollapse. This pause and turnaround creates a caustic dShandarin & Zeldovich , 
19891) at what will be called the second turnaround radius, ri t . 

At yet earlier times, regions of yet higher density collapsed and resulted in caus- 
tics of higher order turnarounds. In the real world, departur es from spherical sym- 
metry would result in violent relaxation (iLvnden-Belll 1 1 9671) . departures from radial 
orbits, and blurring of the caustics. The collapsed region will begin to evolve toward 
energy equipartition. Such regions will be referred to as 'quasi-virialized'. 

These theoretical musings provide a basis for a definition of groups. It will be 
argued that radii of second turnaround and, more problematically of first turnaround, 
are observable. 

The caustic of second turnaround for a group can potentially be defined by two 
features: a density disco ntinuity and a veloci ty dispersion discontinuity. These fea- 
tures were displayed by iBertschingerl (119851) in the simple case of spherical infall 
of an initial top hat density excess. In this simple case, density falls off as a power 
law with radius within r 2t and velocity dispersions are large. Beyond rit, the density 
is expected to take a downward step and velocities locally are expected to manifest 
coherent infall with modest dispersion. The same salient feature s are recovered with 



N-bo dy simulations based on the same spherical infall model ( IVogelsbereer et al 



20091) .' These features might be hard to identify in observational information. The 
density step might be obscured by noise from line-of-sight contamination. The step 
from large to small velocity dispersions can easily be confused by the limitation 
of access to only radial velocities and lack of precise distances that would distin- 
guish backside from frontside infall. The situation with many groups will be messy 
but there are relatively clean cases where the desired information can be recovered. 
For example one can identify the caustics in redshift-radius plots for rich clusters 



(IDiaferio & Gellen.ll997t iRines et all 120031) . The challenge given attention here is 



to find analogous features in groups with much smaller populations. The challenge 
can be met if the entities are sufficiently isolated or so nearby that their three dimen- 
sional structure can be discerned. 

The musings suggest that there should be interesting scaling relations between 
groups. In the spherical approximation, the density at r 2t is the same for all groups 



M(r 2t ) 



constant. 



(5) 



With the approximation M(r) « r 
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M(r 2t )=M v (r 2t /r g ) (6) 

so 

— • — = = constant. (7) 

Hence 

r 2 ,oc (7 p ocMv /3 . (8) 

In subsequent sections, data will be presented that confirm these relationships and 
provide correlation zero points. 

There is the prospect that equivalent scaling relationships can be found based on 
the first turnaround radius, r\ t . The interest here is that the ratio r^jr^t is depen- 
dent on the dark energy content of the universe. Collapse conditions at a given time 
depend almost exclusively on the matter density but dark energy affects the clock, 
hence the timing between collapse events. In a flat universe with 70% of the cosmo- 
logical density in vacuum energy then r\ t jri t = 3.3 whereas in a flat universe with 
the critical density in matter then r\,/r2 t = 3.7 (Gary Mamon, private communica- 
tion). 

Finally it is to be entertained if there are observational manifestations of the 
'zero gravity' surface at a radius r zg around a halo. Only material wit hin this surface 



will e ver participate in collapse onto the halo. It has been suggested ( IChernin et al 



20091) that the region of the zero gravity surface is expected to be evacuated and 



can be identified by this property. The importance of this feature depends on the 
properties of dark energy. 



2 From Big to Small 

The questions to be addressed are whether there are observational features that dis- 
tinguish the quasi-virialized and infall regions of halos and, if so, how these features 
scale with halo mass. Most of the discussion in this paper will focus on modest to 
very small groups of galaxies, the visible manifestations for the overwhelming ma- 
jority of identifiable halos. The sites of intermediate mass halos, those in the range 
5 • 10 12 — 10 14 M Q , have been investigated with a wide field imaging survey at the 
Canada-France-Hawaii Telescope (CFHT), initially using the 0.3 sq. deg . 12K CCD 



camera and then the 1 sq. deg. Megacam detector. ([Jacobs et aU 12009b . Follow up 
spectroscopy was undertaken with Subaru and Keck telescopes. The study is ex- 
tended to very small halos in the range 10 11 - 5- 1O 12 M by giving consideration 
to the region within 4 Mpc where, outside the zone of obscu ration, almost every 



galax y brighter than Mb = —11 has probably been identified ( Karachentsev et al 



20041) . Almost all of these nearby galaxies have now been sufficiently observed with 
Hubble Space Telescope (HST) that an accurate distance is available from a mea- 
surement of the luminosity of the tip of the red giant branch, the TRGB method 
There is now detailed information on the grouping properties of galaxies down to 
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Fig. 1 Radial distributions of galaxies in the Coma and Virgo clusters. The density discontinuities 
at 3 Mpc and 2.0 Mpc respectively are associated with the caustics of second turnaround. 



the scales of associations of dwarfs (ITullv et al. , 2006). Page limitations prevent a 
detailed discussion of the many individual group environments that have been given 
attention. The following is an expurgated version of a study that will be flushed out 
in a journal article. 

The Coma Cluster rests in the nearest massive and evolved halo at 2 • 10 15 M and 
the Virgo Cluster is the dominant structure in the Local Supercluster at 8 • 10 14 M . 
Each is continuing to grow by infall along filaments. The dimensions orthogonal to 
these filaments are well defined. It is shown in Figure[T]that density discontinuities 
give good definition of the caustics of second turnaround. In the case of the Virgo 
Cluster where there is improving distance resolution the location of the caustic is 
supported by velocity information. 

The CFHT wide field imaging survey gave attention to thr ee groups with masses 



in the range 2 - 8 ■ 1 u M P) at distances 25 - 29 Mpc dMahdavi et all 12005 



Trentham et al. . 120061: ITullv & Trenthami 12008). The evidence for the detection of 



second turnaround caustics is shown for two of these cases in Figure [2] Two of 
these intermediate mass halos are dynamically evolved, with central dominant el- 
lipticals and gas-poor companions, while the third presents a situation with mixed 
morphologies. 

In an effort to understand what is going on in low density environments, a sub- 
stantial area around the g alaxy NGC 1023 was imaged with the CFHT cameras 
dTrentham & TullvL 120091) . The central halo has a mass of 6- 10 12 M Q . The detec- 
tion of the second turnaround caustic and the morphological segregation of galaxies 
is seen in Figure[3] 

The very complete information on the three dimensional distribution of galaxies 
within 4 Mpc lets us look at very small halos. The largest of the nearby halos is 
centered on Centaurus A and has a mass of 1 • 10 13 Mq. It is seen in Figure [4] that 
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Fig. 2 Radial distributions of galaxies in the NGC 5846 and N5353/4 groups. The density discon- 
tinuities at 840 kpc and 580 kpc respectively are associated with the caustics of second turnaround. 
The velocity discontinuity at the caustic is shown in the case of NGC 5846. 
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Fig. 3 Left: Galaxies in the vicinity of NGC 1023. Red: early; blue: late; green: transition. Irreg- 
ular box: region of CFHT wide field imaging survey. Dotted circles: inferred caustics of second 
turnaround for the 3 dominant halos. Right: Radial distribution of galaxies with respect to NGC 
1023. The caustic of second turnaround is associated with the discontinuity at 350 kpc. 
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Fig. 4 Left: Distribution of galaxies in the region of Cen A and M83. The slanted dotted line 
is a projection of Galactic latitude +15°. Accurate distances reveal that Cen A and M83 are the 
dominant galaxies in two distinct groups. Right: The cumulative number of early (solid red) and 
late (dotted blue) galaxies as a function of projected radius from Cen A. 



the Cen A and M83 halos are discrete and the Cen A halo is host to a substantial 
family of gas-poor dwarfs (IKarachentsev et al. ■ I2002bll2007h . 

Figure [5] gives attention to the historical Local Group, an environment with two 
discrete halos with masses 1 - 2 ■ 10 12 M Q within a common infall zone. There is 
clear morphological segregation between the regimes inside and outside the second 
turnaround caustics. There is reasonable definition of the surface of first turnaround 



(IKarachentsev et aUl2009l) 



The M8 1 Group provides a relatively clean situation (IKarachentsev et aUl2002al) . 



The central halo of 4 • 10 M Q is characterized by a large dispersion of velocities 
and it is surrounded by a pattern of infall extending to a first turnaround or zero 
velocity surface at 1.3 kpc. See Figure [6] The second panel of this figure illustrates 
the more uncertain situation with the ~ 1 ■ 1 12 M© halo centered on NGC 253 
dJerien et al. . 1 1998t IKarachentsev et all 12003b . 

Finally, Figure [7] illustrates the velocity patterns around two halos in the range 
4 — 7 • 10 11 Mm that on ly contain minor galaxies with luminosities of the LMC or 
less ( Tullv et al l l2006h . These are interesting because luminosities are so small that 
mass to light ratios are high. It is suspected that smaller halos become unidentifiable 
because of a lack of galaxy signposts. 
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Fig. 5 Left: Two projections of the distribution of galaxies around the Local Group. Blue: late 
types; red: early types; green: transition. The slanted dotted lines are projections of Galactic lati- 
tude ±15°. The inner dotted circles approximate the radii of second turnaround for the M31 and 
Milky Way halos. The outer heavy dashed circle approximates the first turnaround radius for the 
combined halos. A light dashed circle encloses the separate NGC 3109 association. Right: Veloci- 
ties as a function of distance from the nearer of M3 1 or the Milky Way. Inside the caustic of second 
turnaround at roughly 270 kpc around each of the major galaxies the velocity dispersion is large 
and galaxies generally have early types. Outside second turnaround most galaxies are late types 
and there is a pattern of infall with the first turnaround located at 940 kpc. 



3 Scaling Relations 

It was anticipated in the introduction that the two independent variables r2 t and a p 
would be directly correlated. The relationship is seen in the left panel of Figure [8] 
The constant of proportionality is found empirically to be 350 km s _1 Mpc . 

Two interesting relations are shown in the right panel of Figure [8] and the left 
panel of Figure|9] Each involves the virial mass calculated according to 

M v = -^i (9) 
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Fig. 6 Infall patterns onto M81 (Left) and NGC 253 (Right). With M81 the differentiation between 
the region of high velocity dispersion within the second turnaround and the region of infall is clear. 
The situation is less clear in the case of the smaller halo containing NGC 253. 
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Fig. 7 Velocities for galaxies in two associations of dwarfs with halo masses below I0 :: M . . 



where the virial radius r g was defined in Eq.Q]so is calculated independently from 

rit- What is being plotted at the right in Fig.[H]is r^t versus Gp^r^ 3 . 
The derived correlation is 

rn = 0.200 Mj 1 ^ 3 Mpc (10) 

if the virial mass is in units of 10 12 M . The tight correlation implies a close cor- 
relation between r2t and r g . It follows that the correlation seen at the left in Fig. [9] 
would exist. In this plot the slope is fixed and the zero-point is determined by the 
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Fig. 8 Left: Correlation between projected second turnaround radius and line-of-sight velocity 
dispersion. Right: Correlation between projected second turnaround radius and halo virial mass. 
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Fig. 9 Left: Correlation between the halo virial mass and line-of-sight velocity dispersion. Right: 
Number of galaxies as a function of the mass of the parent halo, separated between dwarfs (red 
points and solid line fit) and giants (grey points and dashed fit). The separation between giants and 
dwarfs is made at Mr = — 17. 



relations seen in the two panels of Fig. [8] The mass-dispersion relation for groups 
obeys the law 

M v = 2.5 x 10 6 (11) 



where mass M v is in solar units and velocity dispersion a p is in km s~ 
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4 The Radius of First Turnaround 



The best opportunities to measure r\ t , the radius of first turnaround, are presented 
by the groups w ithin 4 Mpc including the Local Gro up. An initial discussion by 



Sandage dl986h was carried forward most recently by iKarachentsev et al.l d2009l) 



The situation in the Local Group is shown in Figure [5] The value of r\ t = 940 kpc 
enclosing a mass of 2 x 10 12 M Q seems reasonably constrained to an uncertainty 
of 10%. However spherical symmetry is a poor approximation for the zero velocity 
surface of the Local Group because of the dumbbell distribution of mass between 
M31 and the Milky Way. Moreover, only 5 galaxies are well placed to define the 
zero velocity surface. 

A second relatively clean opportunity to measure the surface of first turnaround 
is presented by the pattern of infall around M81 seen in Fig. [6] The zero velocity 
surface lies at 1.3 Mpc from the center of a 4 x 10 12 M© halo. 

This limited information nonetheless gives us a preliminary calibration of the 
relationship between r\ t and the mass within this radius: 

n, = 0.78 Mjf Mpc. (12) 

An interesting potential confirmation of the validity of this relation could come 
from the definition of the zero velocity surface around the Virgo Cluster. It is an- 
ticipated to lie at r\, ~ 7 Mpc , a val ue in accordance with long-standing models 



of Virgo infall (T ully & ShayaL 119841) . A dense grid of galaxy distances across the 
Local Supercluster is needed to resolve the issue. 

The relationship between a specific radius like rn or r\ t and the mass internal 
to the specific radius is virtually independent of the dark energy content of the uni- 
verse. However, because of the relative differences with respect to the age of the 
universe between the onset of collapse that is associated with r^t and the current age 
associated with r\ t , the ratio r\ t jri t has a dependence on the dark energy content. 

Currently the best estimate of the ratio of these radii is based on observations 
of the M81 halo. From the data used to generate the left panel of Fig. [6] a tenta- 
tive ratio r\ t /r2t = 3.5 ± 0.4 is found where the error estimate includes random but 
not systematic effects. At present the errors encompass the full range of reasonable 
cosmological models. 

Another quantity of potential interest is the surface of zero gravity, the surface 
around an overdense region that is at the limit of what is destined to ever collapse 
onto the overdensity. This surface is well specified as a function of cosmological 
model in the case of spherical collapse. Wit h £2j\ = 0.7 in a flat universe then rzc ~ 



1.4ri; (IPeirani & de Freitas PachecoL 120081) . However given the close proximity of 



separate collapsed halos, often with more than one halo within an infall region and 
separate infall regions close to overlapping, the assumption of spherical collapse is 
generally a bad approximation. It must be a very bad approximation with rzG- 
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Groups of both the early and late types were studied during the CFHT Megacam 
imaging campaign. The completeness limit for membership studies was typically 
Mr = —11 or fainter. In an effort to study the properties and possible variations 
in the luminosity function, galaxies were separated into two bins: 'giants' with 
Mr < — 17 and 'dwarfs' with — 17 < Mr < — 11. Small but significant variations 
in luminosity functions were found. The variations could be characterized by varia- 
tions in the ratio of dwarfs to giants. The groups of predominantly early types have 
larger dwarf/giant ratios. An interesting result is found if the number of giants or 
dwarfs is plotted against the group virial mass. The correlation is poor with giants 
but very pronounced with dwarfs. The result is shown in Figure|9] The straight line 
describes the relation 

N d =5.2M° 2 9l±Qn (13) 

where N c / is the number of dwarfs with — 17 < —Mr < — 1 1 in the group. The slope 
is consistent within the errors with unity. The number of dwarfs per unit halo mass 
is roughly constant. A count of the number of dwarfs is a good measure of the mass 
of a group halo. 

It is evident from Figs. [3] |U and [5] that the morphological separation with en- 
vironment familiar from studies of rich clusters is also seen in small groups. The 
caustic of second turnaround is a clear marker of the transition between preponder- 
ant morphological types. Secondary galaxies must loose their supplies of cold gas 
within of order the halo crossing time of 

t x = ^J&2i = 2.5 xlO 9 years. (14) 



6 Summary 

Halos can be observationally delineated by a density drop and a transition between 
dispersed orbits and infall at the radii of second turnaround, r 2t - The scaling rela- 
tions that are theoretically anticipated to exist over a wide range of halo masses are 
confirmed from observations of nearby collapsed regions. Specifically: 

• Over the mass range from the halo of M3 1 to the Coma Cluster there is a cor- 
relation between projected second turnaround radius R 2t and line-of-sight velocity 
dispersion a p : o p /R2t = 350 km s _1 Mpc -1 . 

• Over the same mass range, there is the correlation between R 2t and the virial 
mass M12 measured in units of 10 12 M & : R 2t = 0.200(M 12 ) 1//3 Mpc. 

• These two correlations can be combined to give the alternative, but not inde- 
pendent, relation: M V /M Q = 2.5 x I0 6 a p . 

Dwarf galaxies serve to identify collapsed halos because they are numerous. 
Dwarfs and giants alike that are located within the radii of second turnaround tend 
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to be gas-poor 'early' types while galaxies of all sizes outside the second turnaround 
tend to be gas-rich iate' types. The transition between predominantly early and late 
populations is another indicator of the radius of second turnaround. 

Dwarf galaxies delineate collapsed halos in a way that is surprisingly quantita- 
tive. If the number of dwarf galaxies with — 17 < Mr < —11 are counted, then the 
number of these dwarfs in a halo depends linearly on the mass of the halo. Count 
the number of dwarfs and one has a measure of the mass of the halo. 

There are differences in the luminosity function of galaxies with environment 
(Trentham & Tullv, 120091) . More dynamically evolved environments have somewhat 
steeper faint end slopes. However, it appears that the reason is not more dwarfs per 
unit mass in evolved environments but fewer intermediate luminosity systems in the 
vicinity of L*, the luminosity that characterizes the exponential cutoff from the faint 
end power law distribution. It is suspected that intermediate luminosity systems are 
being lost through mergers with the central dominant galaxy. At faint luminosities, 
the production and depletion mechanisms are such that the number of dwarfs per 
unit halo mass remains roughly constant. 
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